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ABSTRACT 

It is proposed that the luminosity function, the comoving-frame spectral correlations and distri- 
butions of cosmological Long-duration Gamma-Ray Bursts (LGRBs) may be very well described as 
multivariate log-normal distribution. This result is based on careful selection, analysis and model- 
ing of the spectral parameters of LGRBs in the largest catalog of Gamma-Ray Bursts available to 
date: 2130 BATSE GRBs, while taking into account the detection threshold and possible selection 
effects on observational data. Constraints on the joint quadru-variate distribution of the isotropic 
peak luminosity (Lj so ), the total isotropic emission (Ei so ), the comoving-frame time- integrated spec- 
tral peak energy (E P:Z ) and the comoving-frame duration (Tqq.z) of LGRBs are derived. Extensive 
goodness-of-fit tests are presented to gauge the consistency of the model predictions with BATSE and 
other gamma-ray instruments data. The presented analysis provides evidence for a relatively large 
fraction of LGRBs that have been missed by BATSE detector with total isotropic emissions (Ei SO ) 
extending down to ~ 10 49 [erg] and observed spectral peak energies (E p ) as low as ~ 5 [ifeV]. The 
intrinsic correlations among these spectral parameters are discussed. Most significantly, the model 
predicts a fairly strong but significant correlation (Pearson's correlation coefficient p — 0.58 ± 0.04, 
at > 14c) between Ei SO & E p>z of LGRBs. Also predicted are strong correlations of Li SO & Ei SO 
with Tgo jZ and a moderate correlation between Li SO & E PjZ of LGRBs. The strength and significance 
of the correlations found, encourage the search for underlying mechanisms, though undermines their 
capabilities as probes of Dark Energy's equation of state at high redshifts. Corroborating recent 
reports on the apparent observed flattening in the Number-Intensity distribution (log(A) — log(P) 
diagram) of LGRBs, the model further extends the reported flattening to a turnover at the bolometric 
1-second peak energy flux P^i ~ 5 x 10~ 8 [erg/s/cm 2 ], corresponding to 1-second peak photon flux 
-P50-300 ~ 0.1 [ph/s/cm 2 ] in BATSE detection energy range: 50 - 300 [KeV]. Although the lack of 
known redshifts for BATSE events limits our ability to set constraints on LGRBs rate, the presented 
analysis favors - but does not necessitate - a cosmic rate for BATSE LGRBs that traces the metallicity 
evolution consistent with a cutoff Z j 1 Zq ~ 0.2 — 0.5, assuming no luminosity-redshift evolution. The 
proposed multivariate model is superior to previously used univariate or multivariate models in the 
GRB literature for its ability to constrain the univariate distributions and quantify the correlations 
of the spectral parameters of LGRBs fully and simultaneously with minimal free parameters. 
Subject headings: gamma rays: bursts — methods: statistical — Gamma-rays: general 



1. INTRODUCTION 

Ever since the discovery of the first G amma- Ray Burst 
(GRB ) by the Vela Satellites in 1967 ( |Klebesadel et al. 
|1973[ ) , there has been tremendous effort and attempts to 
constrain the energetics, luminosity function and the un- 
derlying mechanism re sponsible for these events. E arly 
observations of Konus dMazets fc Golenetskii 1981) and 



Ginga (Fenimore et al. 1988; Nishimura 1988) gamma/X- 
ray instruments suggested a possible link between GRBs 
and Neutron Stars with output energy ranges of the 
order of ~ 10 40 [erg]. With the launch of Compton 
Gamma-Ray Observatory (CGRO), the Burst & Tran- 
sient Alert Telescope (BATSE) onboard CGRO dramat- 
ically changed the u nderstanding of GRB s. While pre- 
vious catalogs (e.g., 



1987) indicated an 



Atteia et al 

isotropic distribution of GRB sources, the BATSE ob 
servations extended this isotropy down to the weakest 
bursts. The non-homogenous (e.g.jFenimore et al.|1993|) 
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and i sotropic spati al event distributio n (e.g ., IMeegar 
et al.||1992| |Briggs 1993] |Fishman et al.|[l994| provided 
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for the first time, strong support for a cosmological ver- 
sus galactic origin of GRBs, undermining neutron stars 
in the local universe as the potential candidates for some 
- if not all - classes of gamma-ray events. Furthermore, 
the joint duration-hardness distribution of GRBs pro- 
vided a direct evidence for at least two separate classes 
of Gamma-Ray Bursts: Long-soft vs. short-hard (e.g., 
Kouveliotou et al.]|1993 also Figure [I] here). 



The possibility of a cosmologicaU origin for GRBs, 
indicated an enor mous output energy on the order of 
~ 10 51 [erg] (e.g., | Dermer||l99"2| ). Nevertheless, an accu- 
rate description or the GRB Luminosity Function (LF) 
also required a knowledge of GRB cosmic rate, an infor- 
mation that could not be extracted from BATSE observa- 
tions alone. This became possible only with the launch of 
the It alian-Dutch X-ray satellite BeppoSax ( Boella et al. 



19971 and the identification of t he first GRB with fi rmTy 



measured cosmological redshift ( jMetzger et al. 1997) that 
marked the beginning of the afterglow era in the field 
of Gamma-Ray Burs ts. The launch of SWIFT satellite 
(Gehrels et al. 20041 was another milestone that revo- 
lutionized the study of GRBs by facilitating the X-ray 
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afterglow observations ( Burrows et al.|2005 1 and further 
ground-based follow-ups for redshilt measurement. 

Alongside the observational triumphs over a few 
decades, several theoretical models have stood up against 
the rivals based on the available evidence an d GRB data. 
Most prominentl y, the Collapsar model (e.g. , |MacFadyen| 
fc Woosley|1999 ) has been relatively successful in linking 



the Long-duration class of Gamma-Ray Bursts (LGRBs) 
to the final stages in the lives of massive stars while the 
Short-duration class of bursts (SGRBs) is generally at- 
tribut ed to the c oalescence of compact binary systems 
(e.g., Nakar 2007 and references therein). Further refine- 
ment oFtEiepotential candidates as the progenitors and 
the emission mechanism for both classes, requires more 
rigorous analysis of observational data in all possible en- 
ergy frequencies. In particular, the prompt gamma-ray 
emission of LGRBs has been subject of intense observa- 
tional and theoretical studies: 

Beginning with BATSE observations, numerous au- 
thors have examined the prompt emission of LGRBs 
searching for potential unde rlying correlations among 



the spectral parameters (e.g., |Nemiroff et al.||1994| Fen 
imore et al.||1995} |Mallozzi et al.||1995[ |Petrosian fc Lee 



1996| |Brainerd|1997||Dezalay et al.|19971|P ctrosian et al. 
1999] |Lloyd et al.||2000[ |Norris et al.||2005p . The lack 
of known redshitts for BATSE events and poor knowl- 
edge of LGRB cosmic rates, however, strongly limited 
the prediction power of such analyses. Instead, the first 
direct evidence for potential correlations and constraints 
on the distributions of the prompt emission spectral pa- 
rameters came with a few LGRBs detected by BATSE, 
BeppoSax, IPN or HETE-II satellites with measured red- 



shifts (e.g., Reichart & Lamb 2001; Amati et al. 2002 
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however, have 
handful of events with high Signal-to-Noise Ratio (SNR) 
required for spectral analysis with afterglows sufficiently 
bright for redshift measurement, arguing that the pro- 
posed joint distributions of the spectral properties, do 
not r epresent the entire underlying population of LGRBs 



(e.g.,|Band fc Preece 2005 Nakar fc Piran 20051 |Butler| 
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20111). Responding to criticisms, attempts were made 



to model the effects of different gamma-ray instruments' 



analysis (e.g. 
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Despite significant progress, the difficulties in model- 
ing the complex detector threshold on the multivariate 
distribution of the prompt emission spectral properties 
and the lack of a sufficiently large sample of uniformly 
detected LGRBs has led the community to focus on indi- 
vidual spectral va riables, most importantly the luminos- 
ity function (e.g ., Petrosian||1993| |Schmidt|[l9 99; Kom- 
niers et al.|2000[|Kumar fc Piran|2000[|Porciani fc Madau 



2001 



Sethi & Bhargavi 2001; Schmidt 20011 |Stern"et~aT 



20021 ITIuetta et al.|[2005| |Salvaterra fc Ghincarini||2007 
S alva terra et al 
Wanderman & 



20091 |Schmidt|20091 ICampisi et al 12010 



nanpOlO; Salvaterra et al. 2012 file 



cently, |Butler et al.| ( |2010[ hereafter BlOj presented an 
elaborate multivariate analysis of SWIFT LGRBs, in- 



cluding the potential correlations among three spectral 
variables: total isotropic energy emission (Ei SO [erg]), a 
definition of duration (T r ^§ [s]), and the spectral peak 
energy (E PtZ [KeV]) of the bursts. Focusing their anal- 
ysis on the E iso — E pz — T r 4 5 interrelation, B10 find a 
strong and significant correlation - but with a broad scat- 
ter - between the isotropic emission and the peak energy 
of the SWIFT LGRBs. Also realized by B10, is the pos- 
sibility of a positive - and perhaps strong - correlation 
between the duration and the the total energy output 
of the bursts. Moreover, to accommodate the potential 
existence of a large population of sub-luminous events, 
a broken power-law luminosity function (LF) is used by 
B10. 

It is known that the energy budget of Gamma-Ray 
Bursts must be limited and a turnover in the LF at low- 
luminosity tail of the population is expected. However, 
the choice of the broken power-law as the candidate LF 
is justified by the fact that current observational data 
cannot constrain this turnover point - expected to be 
far below the detection threshold of current gamma-ray 
instruments. The existence of a turnover point in the 
LF can have important clues for the underlying physics 
of LGRBs. The low-luminosity tail of the LF in such 
case would be the result of the convolution of the Stellar 
Mass distribution with the LGRB rate as a function of 
the properties of massive dying stars, imposed by the 
mechanism. 

Motivated by the search for the potential shape of the 
LGRB LF at the low-luminosity and the flurry of recent 
reports on the possible existence of correlations among 
LGRB spectral variables, the authors of this manuscript 
have considered a wide variety of different statistical 
models that could incorporate and explain all observed 
correlations and spectral distributions, while preserving 
the prediction power of the model for the spectral prop- 
erties of potentially large fraction of LGRBs that could 
go undetected by the current gamma-ray detectors. Here 
it will be shown that it is indeed possible to construct 
LGRB world models capable of describing most (if not 
all) LGRB prompt-emission spectral properties observed 
in the current data catalogs. Towards this, the presented 
analysis is mainly focused on the largest catalog of GRBs 
available to date: BATSE catalog of 2130 GRBs. The 
authors show that there is still an enormous amount of 
information buried in BATSE observations that need to 
be dug out by GRB researchers. 



In the following sections, an example of such data min- 
ing on BATSE catalog will be presented: Section 2.1 
presents an elaborate method for classifying Gamma- Ray 
Bursts into the two known subclasses of S hort and Long 
duration (SGRBs & LGRBs). In Section |2~2| an LGRB 
world model capable of describing BATSE observations 
is presented, followed by a discussion of t he procedure for 
fitting the model to data in Section [2~3) Univariate and 
multivariate goodness-of-fit tests are performed to e nsur e 
the model predicts BATSE data accurately (Section 2.4). 
The results from model fitting and implications on the 
multivariate distribution of the prompt emission spectral 
properties are discussed in Section [3] A summary of the 
analysis and important outcomes ofthe proposed LGRB 
world model will be presented in Section [4] 



LGRB World Model 
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2. GRB WORLD MODEL 

2.1. Sample Selection 

Depending on their detection criteria, some gamma- 
ray detectors might facilitate the detection of one class of 
bursts over the others. For example, the specific detector 
sensitivity of the Busrt Alert Telesc ope (BAT) onboard 
SWIFT satellite (|Gehrels et al.|2004[ ) results in better de - 



2003 2006) 



tections of LGRBs over SGRBs (e.g., |Band 
Therefore, a simple GRB classification method, such as 
a cutoff in the observed duration distributions of GRBs 
(Tgo ~ 3) generally results in LGRB sam ples that are 
minim ally contaminated by SGRBs (e.g., 
2010|. Compared to SWIFT, however, B 



Butler et al. 
ATSE Large 

Area Detectors (LAD) had an increased sensitivity to 
SGRBs. Knowing that many spectral properties of Long 
& Short -duration GRBs, most importantly Tgo, over- 
lap, identification of LGRBs in the BATSE catalog solely 
based on their Tgo will likely result in significant number 
of misclassifications (e.g., Figure|Tj center panel). 

Thus, to ensure a correct analysis of long-duration class 
of GRBs, it is first necessary to collect a pure bias-free 
sample of events, all belonging to LGRB category. This 
goal is achieved by testing extensive varieties of classi- 
fication and clustering method s, most important l y, the 
fuzzy clustering algorithms of Rou sseeuw et al l (|1996f 



and t he method of fuzzy C-means (|Dunn||1973J [Bezde 



3) 



1981 1. Each BATSE-catalog GRB is assigned a probabil 
ity (i.e. class coefficient) of belonging to LGRB (versus 
SGRB) population according to the choice of the clus- 
tering algorithm and the set of GRB spectral variables 
used. This can be any combinations of the peak flux 
(-P50-300 [photon/ s/ cm 2 ]) in BATSE detection energy 
range, in three different time scales: 64, 256 & 1024 [ms]; 
bolometric fluence (Sboi \erg/c m 2 ]); the spectral peak 
energy (E p [KeV]) estimates by |Shahmoradi fc Nemiroffj 
(2010); duration (Tg [s]) and the Fluence-to-Feak-Fiux 
Ratio (FPR [s]). Then GRBs with LGRB class coef- 
ficient > 0.5 are flagged as Long-duration class bursts. 
Overall, the fuzzy C-means classification method on the 
two GRB variables E p & T go is preferred over other clus- 
tering methods and sets of GRB variables (c.f. Appendix 
[A| . This leads to the selection of 1376 events as LGRBs 
out of 1966 BATSE GRBs having measured spectral pa- 
rameters mentioned above. El 

As a further safety check to ensure minimal contami- 
nation of the sample by SGRBs, the lightcurves of 291 
bursts among 1966 BATSE GRBs with LGRB class coef- 
ficients in the range of 0.3 — 0.7 are visually inspected in 
the four main energy channels of BATSE Large Area De- 
tectors (LAD). This leads to reclassification of 17 events 
(originally flagged as LGRB by the clustering algorithm) 
to potentially SGRB or Soft Gamma Repeater (SGR) 
events, and reclassification of 7 events (originally flagged 
as SGRB by the clustering algorithm) to LGRBs. The 
result is a reduction in the size of the original LGRB 
sample from 1376 to 1366 (Table [2]). 

1 Data for 1966 BATSE GRBs with firmly measured peak 
flux, fluence and duration are taken from The BATSE Gamma 
Ray Burst Catalogs: http://www.batse.msfc.nasa.gov/batse/ 
grb/catalog/ The spectral peak energy (E„) esti mates of these 
events are taken from Shahmoradi & Nemiroff ( 2010 ), also available 
for download at: https://sites.google.com/site/amshportal/ 
research/aca/ in- the-news/lgrb-wor ld-model 
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Fig. 1.— Classification of 1966 BATSE LGRBs according to 
the mos t sui table clustering algorithm and set of GRB variables 
(Section [2Tj: Fuzzy C -means classification on E p &c Tgo- Red & 
blue colors represent LGRB & SGRB classes respectively, in all 
three plots. Top: The joint Tgo — E p distribution. The uncer- 
tainties in LGRBs are derived from the Empirical Bayes model 
discussed in Appendix [C] Center. Tgo di stribution. Bottom: E v 
distrib ution. E p estimates are taken from [Shahmoradi & N emiroff 
J2010[l. Compare t his plot to the plot of Figure (13 ) of Shahmoradi 
|fc INem iroff (2010), where the entire univariate E p distribution ot 
rJAlHt; liltiss were fit by a 2-component Gaussian mixture. 



4 



Shahmoradi 



2.2. Model Construction 

The goal of the presented analysis is to derive a mul- 
tivariate model that is capable of reproducing the obser- 
vational data of 1366 BATSE LGRBs. Examples of mul- 
tivariate treatment of LGRB data are rare in Gamma- 
Ray Burst literature, with the most recent (and perhaps 
the only) example of such work presented by BIO. Con- 
versely, many authors have focused primarily on the uni- 
variate distribution of the spectral parameters, most im- 
portantly on the luminosity function (LF). A variety of 
univariate models have been proposed as the LGRB LF 
and fit to data by approxima ting the compl ex detector 
threshold as a step function ( |Schmidt 19991 or an effi- 
ciency grid (e . g. the four-interval efficiency modeling of 
|Guetta et al.| ( |2005| )) or by other approximation meth- 
ods. A more accurate modeling of the LF, however, re- 
quires at least two LGRB observable incorporated in the 
model: the bolometric peak flux (Pboi) and the observed 
peak energy (E p ). The parameter E p is required, since 
most gamma-ray detectors are photon counters, a quan- 
tity that depends on not only P^oi but also E p of the 
burst. This leads to the requirement of using a bivariate 
distribution as the minimum acceptable model to begin 
with, for the purpose of constraining the LF. The choice 
of model can be almost anything, since the current theo- 
ries of LGRBs prompt emission do not set strong limits 
on the shape and range of the luminosity function (or 
any other LGRB spectral variables). 

Here, the multivariate log-normal distribution is pro- 
posed as the simplest natural candidate model capable 
of describing data. The motivation behind this choice 
of model comes from the available observational data 
that closely resembles a joint multivariate log-normal dis- 
tribution for four most widely studied spectral parame- 
ters of LGRBs in the observer-frame: Pboi, Sboi (bolo- 
metric fluence), E p , Tqq: Since most LGRBs originate 
from moderate redshifts z ~ 1 — 3 (a fact known thanks 
to Swift satellite), the convolution of these observer- 
frame parameters with the redshift distribution results 
in negligible variation in the shape of the comoving- 
frame joint distribution of the same LGRB parameters. 
Therefore, the redshift-convoluted 4-Dimensional (4D) 
observer-frame distribution can be well approximated as 
a linear translation from the observer-frame parameter 
space to the comoving-frame parameter space, keeping 
the shape of the distribution almost intact. This im- 
plies that the joint distribution of the intrinsic LGRB 
variables: the isotropic peak luminosity (Li SO ), the to- 
tal isotropic emission (Ei so ), the comoving-frame time- 
integrated spectral peak energy {E PyZ ) and the comoving- 
frame duration (Tg z ) might be indeed well described by 
a multivariate log-normal distribution. 

In general, models with higher nonzero moments than 
the log-normal model can also be considered fo r fitting, 
such as a multivariate Skew-lognormal (e.g., Azzalini 
|1985| . This, however, requires fitting for higher number 
of tree parameters which is practically impossible given 
BATSE data with no available redshift information. 

Following the discussion above, the process of LGRB 
observation can be therefore considered as a non- 
homogeneous Poisson process whose mean rate param- 
eter - the cosmic LGRB differential rate, lZ cosm ic ~ is 
the product of the differential comoving LGRB rate den- 



sity £(z) with ap = 4 dimensional log-normal Probability 
Density Function (pdf), CM , of four LGRB spectral vari- 

with location vector p, 



ables: L; 



Ei sn , E„ ,_ and To, 



and the scale (i.e. covariance) matrix E, 



dN 



dLiso dEi SO dE PiZ dTgo :Z dz 
cx CAf(^Liso, Eis 0l Ep,z> Tqo,z\ & £ 



(1) 



(1+2) ' 

where the factor (1 + z) in the denominator accounts 
for cosmological time dilation and the comoving volume 
element per unit redshift, dv /dz, is, 



dV _ C_ 

"cL7 ~ ~~ 



4nD L 2 (z) 



{1 + zf 



1/2 ' 



(2) 



with Dl being the luminosity distance, 



D L (z) 



C_ 
H~ 



(1- 



dz' 



(i + z') 3 n M + n A 



-1/2 



(3) 



assuming a flat ACDM cosmology, with parameters 
set to h = 0.70, Q, M = 0.27 and fi A = 0.73 (Jarosik et 
al. 2011). Here, C & H = lOO/i [Km/s/MPc] stand for 
the speed of light and the Hubble constant respectively. 

The 4-dimensional log-normal distribution of Equa- 
tion (JlJ , CAf , has an intimate connection to multivariate 
Gaussian distribution in the logarithmic space of LGRB 
observable parameters (c.f. Appendix O). 

One could generalize the LGRB rate of Equation (HI to 
incorporate a redshift evolution of the LGRB variables in 
the form of fJ*i(z) — /io,i+c^ log(l+z), (i = 1, .., 4), where 
a has to be constrained by observational data. This is, 
however, impractical for BATSE data due to unknown 
redshifts, as the fitting results in degenerate values for 
a. Nevertheless, the multivariate analysis of SWIFT 
LGRBs presented by B10 strongly rejects the possibil- 
ity of redshift evolution of the LF, a fact that further 
legitimizes the absence of redshift-luminosity evolution 
in Equation 0. 

As for the the comoving rate density C,{z), it is assumed 
that LGRBs trace the Star Formation Rat e (SFR) in 
the form of a p iecewise power-law function of Hopkins & 
Beacom| (p006| (hereafter HB06), 



' (1 + z)i° 



\{l + z)^ 



z < z 

z Q < z < Z\ 

z> Zi, 



(4) 



with parameters (zq, Z\, 70, 71, 72) set to best-fit values 
(0.97, 4.5, 3.4, -0.3, -7.8) of HB06, also to the best values 
(0.993 , 3.8, 3.3, 0.055, -4.46) of an updated SFR fit byp 
(2008). Alternatively, the bias-corrected redshift distri- 
bution of LGRBs derived from SWIFT data (B10) with 
best-fit parameter values (0.97,4.00,3.14,1.36,-2.92) 

can be employed as C(z). This parameter set is consistent 
with an LGRB rate scenario tracing metallicity-corrected 
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SFR with a cutoff ZI Z^ ~ 0.2-0.5 (Figure (10) & Equa- 
tion (8) in BIO; III (|2008[)). The hypothesis of LGRB 



rate evolving witneosrme metallicity is both predicted 
by the Colla psar model of long-dur ation Gamma-Ray 
Bursts (e.g., |Woosley fc He gcr 2006) and supported by 
observations of LGRBs host galaxies (e.g., |Stanek et al" 



2006[ Levesque et al. 



2010a), although the metallicity- 



rate connection and the presence of a sharp metallicity 
cutoff has been challenged by few recen t host galaxy ob- 
servations (e.g., Levesque et al.|2010c | b| and poss ible un- 
known observational biases (e.g., |Levcsquc 2012]). 

The cosmic LGRB rate, lZ cosm i C , in Equation (fTl), al- 
though quantified correctly, does not represent the ob- 
served rate (K obs ) of LGRBs detected by BATSE LADs, 
unless convolved with an accurate model of BATSE trig- 
ger efficiency, n, as a function of the burst redshift and 
rest-frame spectral parameters, discussed in Appendix |B| 

H-obs = Vi-Lisoi 2 , Tgo i z , z) X 7Z cosm i c (5) 

2.3. Model Fitting 

Now, with a statistical model for the observed rate of 
LGRBs at hand (i.e. Equation^, the best-fit parameters 
can be obtained by the methodof Maximum Likelihood. 
This is done by maximizing the Likelihood function of 
the model, given the observational data, using a vari- 
ant of Metropolis-Hastings Markov Chain Monte Carlo 
(MCMC) algorithm discussed in detail in Appendix [Cj 
As mentioned before, the fitting is performed fo r thre e 
redshift distribution scenarios of HB06, B10 fcp]( |2008[ ). 

It is also known that Tgo of LGRBs are potentially sub- 
ject to estimation biases. To ensure that the reported Tgo 
of BATSE LGRBs do not bias the fitting results for the 
rest of the parameters, model fitting was also performed 
by considering only three spectral variables of BATSE 
LGRBs: the bolometric 1-sec peak flux (Pboi), bolomet- 
ric fluence (Sboi), observed peak energy (E p ), excluding 
the duration (Tqq >z ) variable from the model, thus re- 
ducing the dimension of the model by one. Only after 
the fitting was performed, it became clear that the in- 
clusion of T 90 of BATSE LGRBs in fitting does not sig- 
nificantly affect the resulting best-fit parameters of the 
model. Therefore, here only results from the full model 
fitting are presented, as in Table (l}. 

Due to lack of redshift information for BATSE GRBs, 
the resulting parameters of the model exhibit strong co- 
variations with each other. This is illustrated in the ex- 
ample correlation matrix of LGRB world model in Table 
([3| . All location parameters appear to correlate strongly 
positively with each other, and so do the scale parame- 
ters. The location parameters however negatively corre- 
late with the scale parameters, meaning that an increase 
in the average values of the comoving-frame spectral pa- 
rameters reduces the half- width of the corresponding dis- 
tributions of the variables. In general, it is also observed 
that the correlations among the spectral variables weaken 
with increasing the location parameters. An exception to 
this, is the correlation of Tgo jZ with Lj s0 & Ei SO which 
tends to increase with location parameters. Since an 
excess in the cosmic rates of LGRBs at high redshifts, 
generally results in an increase in the values of location 
parameters, it can be said that 'given BATSE LGRBs 
data, a higher rate for LGRBs at distant universe, gen- 



TABLE 1 

Mean Best-fit parameters of LGRB World Model, 

FOR THREE REDSHIFT DISTRIBUTION SCENARIOS 



Parameter 


HB06 


Li (2008) 




B10 


Redshift Parameters (Equation 


1 




20 


0.97 


0.993 




0.97 


21 


4.5 


3.8 




4.00 


70 


3.4 


3.3 




3.14 


71 


-0.3 


0.0549 




1.36 


72 


-7.8 


-4.46 




-2.92 


Location Parameters 


log(L iso ) 


51.35 ±0.20 


51.50 ±0.19 




51.73 ±0.19 


\og(E lso ) 


51.82 ±0.20 


51.94 ±0.20 




52.03 ±0.21 


log(%,») 


2.43 ±0.05 


2.47 ±0.05 




2.54 ± 0.05 


log(T 90 , z ) 


0.99 ±0.03 


0.96 ± 0.03 




0.80 ± 0.03 


Scale Parameters 


l°g(°"L lso ) 


-0.22 ±0.06 


-0.23 ±0.06 




-0.20 ±0.05 




-0.07 ± 0.03 


-0.07 ± 0.03 




-0.04 ± 0.03 




-0.44 ±0.02 


-0.44 ± 0.02 




-0.44 ± 0.02 




-0.38 ±0.01 


-0.39 ±0.01 




-0.40 ±0.01 


Correlation Coefficients 



PL lso 


-Eiso 


0.93 ±0.01 


0.94 ±0.01 


0.96 ±0.01 




~Ep,z 


0.47 ± 0.07 


0.45 ±0.07 


0.44 ± 0.08 


PL iso 


- T 90,z 


0.52 ±0.08 


0.59 ±0.09 


0.75 ±0.07 


PE i30 


-E P ,, 


0.58 ±0.04 


0.58 ±0.04 


0.59 ±0.04 


PE iso 


- T 90,z 


0.63 ±0.05 


0.66 ± 0.05 


0.74 ± 0.04 




-T90,s 


0.34 ± 0.04 


0.37 ± 0.04 


0.50 ± 0.04 



BATSE LGRB Detection Efficiency (Eqn. A5) 



^thresh 
lg(<Tthresh) 



-0.44 ±0.02 

-0.88 ±0.05 



-0.45 ± 0.02 
-0.90 ± 0.05 



-0.44 ± 0.02 
-0.88 ±0.05 



Note. — The full Markov Chain sampling of the above pa- 
rameters from the 16-dimensional parameter space of the 
likelihood function (Appendix [cl are avail able for down- 
load at https : //sites . google . conT/site/amshportal/research/ 
aca/in-the-news/lgrb-world-model tor each ot the three redshitt 
distributions. 



erally implies weaker Ei so — E p>z & L 
lotions and stronger covariation of Tgo 
other spectral parameters '. 



t-.it Ep^ z corre 
with the three 



2.4. Goodness-of-Fit Tests 

In any statistical fitting problem, perhaps more im- 
portant than the model construction is to provide tests 
showing how good is the model fit to input data. For 
many univariate studies of the GRB LF, this is done by 



1948 lor Pearson's x 



|Kolmogoroff 


1941 


Smirnov 


Fisher 


1924 


tests. In the 



case of multivariate studies (e.g. LJ1UJ, a combination 
of visual inspection of the fitting results, KS test on the 
marginal and bivariate distributions and variants of \ 2 
(e.g. likelihood ratio) tests have been used. 

In general, univariate tests on the marginal distribu- 
tions of multivariate fits provide only necessary - but 
not sufficient - evidence for a good multivariate fit. Al- 
ternatively one could assess the similarity by the use of 
nonparametric multivariate Goodness-of-Fit (GoF) tests. 
Such tests, although exist, have been rarely discussed and 
treated in statistics due to di fficulties in the interpreta- 
tion of the test statistic (e.g., |Peacock||1983 Press et al 



1992 Justel et al.||1997[ ). Ideally, one can always use the 



Pearson's \ l GoF test for any multivariate distribution. 
However, for the special case of BATSE 1366 LGRBs, one 
would need an observed sample consisting of N ^> 1366 



G 
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observations to avoid serious instabiliti es that occur i n 
tests due to small sample sizes (e.g., )Cochran||1954| ). 



X 



To ensure a good fit to the observational data m all 
- and not only univariate - levels of the multivariate 
structure of data, as assessment of similarity can be ob- 
tained by scanning and comparing the model and data 
along their principal axes, in addition to univariate tests 
on the marginal distributions. Although statistically not 
sufficient condition for the multivariate similarity of the 
model prediction to data, this can provide strong evi- 
dence in favor of a good fit, at much higher confidence 
than tests performed only on the marginal distributions. 

Following the lines above, Figure ([2| presents the 
model predictions for marginal distributions of the four 
LGRB spectral variables in the observer-frame. The 
Kolmogorov-Smirnov (KS) test probabilities for the sim- 
ilarity of the model predictions to the marginal distribu- 
tions of BATSE LGRB variables are also reported at the 
top right of each plot. All three redshift-distribution pa- 
rameters of HB06, Li (2008) & BfO (Equation [1 result 
in relatively similar fits to BATSE data in the observer 
frame. Thus, for brevity, only plots for one repre sen- 



Li 2008) are 



tative (median) redshift distribution 
presented. 

As the second level of GoF tests, the joint bivariate 
model predictions are compared to BATSE LGRB data, 
presented in Figures p|) , Q & |5]) . This method of scan- 
ning model and data along the principal axes of the joint 
bivariate distributions can be generalized to trivariate 
and quadruvariate joint distributions. For brevity, how- 
ever, only the bivariate tests are presented. 

3. RESULTS AND DISCUSSION 

It is observed in the plots of Figures S, Q & (|5) 
that the model provides excellent fit to aata, within tnc 
uncertainties caused by random poisson fluctuations in 
the BATSE LGRBs observed rate. These random fluc- 
tuations in BATSE detections are encompassed in each 
graph by the dashed-green lines that represent the 90% 
Confidence Intervals (CI) on BATSE LGRB detections 
(blue solid lines), derived by repeated sampling from the 
model. 

Unfortunately, same methods for a comparison of data 
and model cannot be applied in LGRBs comoving-frame, 
due to lack of redshift information for BATSE sample of 
LGRBs. Nevertheless, a comparison of the model with 
observational data of other instruments, - with measured 
redshifts - can provide clues on the underlying joint dis- 
tribution of LGRBs spectral variables in the co-moving 
frame compared to different instruments' detections, as 
will be done in the following sections. 

3.1. LGRB Luminosity Function & log(A^) — log(P) 
diagram 

The log(A) — log(P) diagram of Gamma-Ray Bursts 
has been subject of numerous studies in the BATSE 
era, primarily for the purpose of finding signatures of 
cosmological- (vs. galactic-) origins in the LGRB rate. 
The cosmological origin of LGRBs is now well estab- 
lished. Nevertheless, the log(A^) — log(P) diagram can 
still provide useful information for future gamma-ray ex- 
periments. 

Figure pi bottom) depicts the prediction of the LGRB 
world model for the traditional log(A) — log(P) diagram 



for 1-second peak photon flux in the BATSE nominal 
detection energy range 50 — 300 [KeV], both for the 
differential (left panel) & the cumulative (right panel) 
LGRB rate. For all three LGRB cosmic rates consid- 
ered in this work - as in Table ([I]) - the differential 
log(A) — log(-P) diagram shows a peak in the rate at 
-P50-300 ~ 0.1[ph/s/cm 2 ]. Such peak in the LGRB rate 
results in a relative flattening at the dim end of the cu- 
mulative log(A) — log(P) diagram, as compared to its 
bright end. This observation has already been reported 
by B10 for SWIFT sample of LGRBs, although an en- 
tirely different luminosity function - a broken power-law 
LF - were used by B10 in their multivariate LGRB world 
model (c.f. Equation (2) in B10). 

On the other hand, the peak in the observer-frame 
LGRB rate translates to three relatively different (at 
~ lex level) peaks in the luminosity function (LF) of 
LGRBs. In general, it is observed that the peak of LF 
(i.e. the average 1-second peak luminosity of LGRBs) 
increases with increasing the cosmic rate of LGRBs at 
high redshift. This effect is well depicted in the Top Left 
panel of Figure j[6| for th e three LGRB cosmic rates con- 
sidered: HB06,O( |2008| , B10. 

Compared to the predictions of BlO's LGRB world 
model (c.f. bottom plot of Figure (6) in B10), the log- 
normal model suggests a lower peak for BATSE LGRBs 
LF (~ 52.3 here, vs. ~ 52.7 in B10) for the same 
redshift distribution of LGRBs. Averaging over the 
three redshift distributions considered, the model pre- 
dicts a dynamic 3<r range of observer-frame brightness 
\og(Pboi[erg/s/cm 2 ]) € [—7.11 ± 2.66] corresponding to 
P bo i[erg/s/cm 2 } € [1.70 x lO" 10 , 3.58 x 10~ 5 ] for LGRBs. 
This translates to an average dynamic 3cr range - in the 
comoving-frame - of log(Li SO [er<7/s]) € [51.53±1.99] cor- 
responding to L iso [erg/s] e [3.46 x 10 49 ,3.38 x 10 53 ]. 

3.2. Isotropic Emission & Peak Energy distributions 

A comparison of the Top Right panel of Figure ([6]) with 
SWIFT observations of LGRBs (e.g. B10, Figure 7 & 8) 
indicates the similarity of the peak rates in BATSE & 
SWIFT samples of total isotropic emission (Ei SO ) distri- 
butions. The distributions for both instruments show 
a similar peak at log(P lso ) - 52.7. The SWIFT E iso 
distribution, however, spans to relatively lower E iso as 
compared to BATSE. This observation is not surprising 
if one takes into account the strong trivariate correla- 
tions of the three LGRB variables: L. LSO , Ei SO , E p , z , all 
of which play role in LGRB detection by SWIFT and 
BATSE. In fact, a comparison of the spectral peak energy 
(E PtZ ) distributions of SWIFT (e.g. Figure (2) in B10) 
with the model's prediction for BATSE LGRBs (Center 
Left panel of Figure [6 ) reveals the relatively high sen- 
sitivity of SWIFT BAT detector to dim soft LGRBs as 
compared to BATSE LADs. Such difference in the two 
detectors has already been discussed frequently by dif- 
ferent authors (e.g. Band 2006). 

In sum, averaging over the three redshift distribu- 
tions considered, the LGRB world model predicts a dy- 
namic 3fi range of observer-frame bolometric LGRB flu- 
ence log(Sboi[erg/cm 2 ]) € [—6.16 ± 3.01] corresponding 
to S bo i[erg/cm 2 } G [6.82 x 10- 10 ,7.01 x 10" 4 ]. This 
translates to an average 3cr range - in the comoving- 
frame - of \og(E iso [erg}) e [51.93 ± 2.71] corresponding 
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Fig. 2. — The marginal distribution predictions (solid blue lines) of the LGRB world model given BATSE detection efficiency, superposed 
on BATSE 1366 LGRB data (red histograms). The solid grey lines represent the model predictions for the entire LGRB population 
(detected and undetected), with no correction for BATSE sky exposure and the beaming factor (/ ). The 90% confidence intervals (dashed 
green lines) represent random poisson fluctuations expected in BATSE LGRB-detection process. Bottom: The differential (right panel) & 
cumulative (left panel) rate of LGRBs as a function of peak photon flux in the BATSE nominal detection energy range (50 — 300 [-R"eV]). 
For a comparison with the SWIFT sample of LGRBs and the proposed multivariate model of B10, the two dashed grey lines - taken from 
Figure (1) in B10 - represent the predictions of the broken power-law luminosity function in Equation (2) of B10, in the observer-frame. The 
Kolmogorov-Smirnov (KS) test probabilities for the goodness-of-fit of the model predictions to BATSE data are reported at the top-right 
corner of each plot. 
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Fig. 3. — Top: BATSE 1366 LGRB data (red dots) superposed on the joint bivariate distribution predictions (black dots) of the LGRB 
world model for BATSE detection efficiency. The background grey dots represent the model predictions for the entire LGRB population 
(detected and undetected). The uncertainties in BATSE LGRBs are derived from the Empirical Bayes model discussed in Appendix [C| 
Center & Bottom: Gauging the goodness-of-fit of the bivariate model to data by scanning and comparing the joint distributions of tnc 
model & BATSE data along their principal axes. The colors bear the same meaning as the top panels, in addition to the 90% confidence 
intervals (dashed green lines) that represent random poisson fluctuations expected in BATSE LGRB-detection process. 
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Fig. 4. — Top: BATSE 1366 LGRB data (red dots) superposed on the joint bivariate distribution predictions (black dots) of the LGRB 
world model for BATSE detection efficiency. The background grey dots represent the model predictions for the entire LGRB pop ulation 
(detected and undetected). The uncertainties in LGRBs are derived from the Empirical Bayes model discussed in Appendix [C] Center 
& Bottom: Gauging the goodness-of-fit of the bivariate model to data by scanning and comparing the joint distributions of the model & 
BATSE data along their principal axes. The colors bear the same meaning as the top panels, in addition to the 90% confidence intervals 
(dashed green lines) that represent random poisson fluctuations expected in BATSE LGRB-detection process. 
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Fig. 5. — Top: BATSE 1366 LGRB data (red dots) superposed on the joint bivariate distribution predictions (black dots) of the LGRB 
world model for BATSE detection efficiency. The background grey dots represent the model predictions for the entire LGRB pop ulation 
(detected and undetected). The uncertainties in LGRBs are derived from the Empirical Bayes model discussed in Appendix [C] Center 
& Bottom: Gauging the goodness-of-fit of the bivariate model to data by scanning and comparing the joint distributions of the model & 
BATSE data along their principal axes. The colors bear the same meaning as the top panels, in addition to the 90% confidence intervals 
(dashed green lines) that represent random poisson fluctuations expected in BATSE LGRB-detection process. 
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Fig. 6. — Top & Center: The marginal distribution predictions (blue dottcd/dashcd/solid lines) of the LGRB world model for BATSE 
LGRBs in the burst comoving frame. The grey dotted/dashed/solid lines — corresponding to the three LGRB redshift distributions 
HB06/Li(2008)/B10 - represent model predictions for the entire LGRB population (detected and undetected), with no correction for 
BATSE sky exposure and the beaming factor (/;,). For a comparison with SWIFT sample of LGRBs, refer to Figures (2), (6) & (7) of B10. 
Bottom: The three differential (right panel) & cumulative (left panel) rates of LGRBs at a given redshift (z). 
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Fig. 7. — The joint bivariato distribution predictions (black dots) 
of the LGRB world model for BATSE detection efficiency, assum- 
ing LGRB rate tracing the cosmic Star Formation Rate (SFR), 
updated by Li (2008). The background grey dots represent the 
model predictions for the entire LGRB population (detected and 
undetected). Superposed on the model predictions are the well- 



known proposed Yonetoku (L;, 



E P ,z) & Amati (£;.. 



.) 



relations. The red lines in each plot represent the best-fit power- 
law relations derived by the corresponding authors. Neither the 
Yonetoku nor the Amati relations appear to be consistent with the 
predicted joint distributions for 1366 BATSE LGRBs (black dots) 
and the entire LGRB population (background grey dots). This 
implies strong bias and selection effects in the detection process, 
redshift determination and spectral analysis of the LGRB samples 
that were used to construct the two relations. 

to E lso [erg] € [1.66 x 10 49 ,4.46 x 10 54 ]. 

As for the spectral peak energy (E p & E P:Z ) distri- 
butions, the model predicts a 3a range of observer-frame 
LGRB spectral peak energy log(E p [KeV}) € [1.93±1.22] 
corresponding to E p [KeV] £ [5,1427]. This trans- 
lates to an average 3rj range - in the comoving-frame 
- of \og(E p z [KeV]) € [2.48 ± 1.12] corresponding to 
E PtZ [KeV] € [23,4006]. 

3.3. Duration distribution 

GRBs are traditionally flagged as long-duration class 
of bursts if their observed durations (Tg ) exceed 2 sec- 
onds (e.g., Kouveliotou et al. 19931. Such classifica- 
tion, however, has been long known to be ambiguous 
close to the cutoff set at T go = 2 [s]. It will be there- 



fore useful to explore how accurate such classification is 
for the entire LGRB population (including non-triggered 
LGRBs). Figure |2] Center Right plot), depicts the un- 
derlying population vs. 1366 BATSE LGRBs observed 
durations (Tgo). As implied by the model, the shape 
of the T 90 distribution of LGRBs is not significantly af- 
fected by the triggering process of BATSE, since both 
BATSE and entire LGRB population distributions show 
a similar peak at Tgo ~ 30 [s] . There is however a slight 
difference (< 0.1 dex) in the predicted observer-frame 
peak of LGRBs T 90 distribution, depending on the un- 
derlying LGRB redshift distribution assumed. The dif- 
ference is magnified to ~ 0.2 dex in the comoving-frame 
(^90,2) duration distribution of LGRBs, for the two ex- 
treme cases of HB06 & B10 redshift distributions. In 
general, a higher LGRB rate at high redshifts (as in the 
case of B10 redshift distribution) results in a shift to 
shorter durations in the duration distribution of LGRBs, 
in both the observer and the comoving frames (Figure pi 
Center Right plot). 

Overall, the model predicts an average dynamic 3tj 
range of observer-frame log(T 9 o[s]) <E [1.47 ± 1.32] cor- 
responding to T 90 [s] e [1.4,620] for LGRBs. This trans- 
lates to an average dynamic 3<7 range of comoving-frame 
log(T 90z [s]) € [0.92 ± 1.24] corresponding to T go z [s] € 
[0.47,145]. 

3.4. Spectral Correlations 

Ever since the launch of SWIFT Gamma-Ray detec- 
tor satellite, the has been a flurry of reports on the dis- 
covery of strong and significant correlations among the 
spectral parameters of LGRBs, most prominently, among 
the comoving-frame spectral peak energy and th e total 
isotropic emission or luminosity of L GRBs (e.g., Amati 
et al.||2008| |Ghirlanda et al 1120081) Despite the lack 
ot measured redshifts tor BATSE GRBs, signatures of 
such correlations had been found by earlier works in the 
BATSE era through careful analysi s of observer-fram e 
spectral properties of LGRBs (e.g., iLloyd et al. 2000). 



Nevertheless, the strength and significance ot these cor- 
relations were undermined by analyzing larger samples of 
BATSE catalog of GRBs (e.g.JNakar & Pir an 2005; JBand| 



& Preece 2005 
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SWIFT sample of GRBs (e.g 



JNcmiroff 
Butler et 



2009,120111) or 
aT pD07f 



2009 



20101, all arguing that the sample ot bursts used to con 



struct the claimed spectral relations is representative of 
only bright-soft LGRBs. These arguments have been 
responded by ot hers (e.g., [G hirlanda et al. 2005, j2008| 
Nava et al.||2008[ c.f. JShahmoradi fc Nemiroflj ( |20ll[ ) for 
a complete history of the debate). 

A multivariate analysis of the BATSE LGRBs data 
that carefully eliminates the potential biases at the de- 
tection process, can therefore greatly help to understand 
the strength and significance of the reported correla- 
tions. Figure [7] shows the predictions of the LGRB world 
model for the two widely discussed spectral correlations: 
Li SO — E p>z (the Yonetoku) and Ei SO — E p>z (the Amati) 
relations. As indicated by the model, a large fraction of 
BATSE LGRBs (and much larger fraction of the entire 
LGRB population) on the dim-hard regions of the plots 
appear to be under-represented by the sample of LGRBs 
used for the construction of these relations. Neverthe- 
less, given the three redshift distributions considered for 
the model, a relatively strong (Pearson's correlation co- 
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efficient pE iso ~E , = 0.58 ± 0.04) and highly significant 
(> 14a) correlation is predicted between Ei SO k E P:Z of 
Long-duration class of GRBs. The slope of the two re- 
lations suggested by the model also diff er significantly 
from the original repo rts of the relations ( Sch aefer||2007| 
Ghirlanda et al. 2008). Generally, in a regression mod- 
eling, it is assumed that there are two independent k 
depenedent variables. In the case of the proposed rela- 
tions, none of the variables is known to depend theoret- 
ically on the other. Therefore, due to the large statisti- 
cally unexplained variances of the two LGRB variables, 
different regression methods, such as Ordinary Least 
Squares: OLS(F|A) k OLS(X|F), result in entirely dif- 
ferent slopes for the relations. The best-fit power-law 
relations and the conditional variances of the regressand 
given the regressor can be easily obtained from the pa- 
rameters of the multiva riate log-normal mo del in Table 
([I]) (e.g., Section 2.11 in [Kutner et al. 

As observed in Table JT]), the mode 
correlation among all 



2004) 



predicts positive 
jGRB spectral variables. The 
isotropic peak luminosity (Li SO ) and the total isotropic 
emission (Ei SO ) appear to be strong indicators of each 
other reciprocally. Surprisingly, it is also observed that 
the comoving-frame duration Tgo,z of LGRBs strongly 
correlates with both Ei SO k and Li SO , The existence of a 
possible positive correlation between the isotropic emis- 
sion and the duration of LGRBs has been implied by the 
analysis of SWIFT LGRBs (B10), though only weakly 
present in there. Such correlations can be enlightening 
for the early studies of time dilation signatures in BATSE 
GRBs (Shahmoradi k Nemiroff, in preparation). A pos- 
itive duration-brightness correlation is also opposite to 
- but not necessarily in contradiction with - the nega- 
tive duration-br ightness correlation in pulses of individ- 



ual GRBs (e.g., Fenimore et al. 1995). Combination of 



the two correlations implies that the number of pulses 
in individual LGRBs should be positively correlated with 
the peak luminosity ( or equivalently, the total isotropic 
emission) of the bursts. This is indeed in qualitative 
agreement with the observed inequality relation between 
the isotropic peak luminosity and the number of pulses 
in SWIFT LGRBs (Figure (6) in |Schaefer|[2007| . The 
strength of the correlations found, encourages the search 
for the underlying physical mechanism that could give 
rise to these relations. This is however, beyond the scope 
of this manuscript (c.f. Rees k M szros 20 05} |Ryde et al. 
20061 [Thompson et al.|2007||Giannios|201ii||Dac[o k Par 



2012 for example discussions} 

It is also worth mentioning that the Tgo duration of 
BATSE LGRBs strongly correlates with the bolometric 
Fluence to bolometric 1-second Peak flux Ratios (FPR), 
with Pearson's correlation coefficient /?fpr-T90 ~ 0.67. 
A comparison of BATSE data with the predictions of 
the model for the bivariate distribution of FPR and Tgo 
is given in Figure [8] Interestingly, the model predicts 
the same correlation strength of /Ofpr-T90 ~ 0.67 for 
the entire LGRB populations, implying that the detec- 
tion process does not bias the FPR — Tgo relation. Such 
strong correlation indicates an underlying intrinsic inter- 
relation between the three spectral variables: P bot , Sboi 
k Tgo, also among their corresponding comoving-frame 
counterparts. In fact, B10 use a variant of this trivariate 
correlation to define an effective peak flux in terms of flu- 
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Fig. 8. — The joint bivariate distribution of duration (Tgo) vs. 
the Ratio of bolometric Fluence to 1-second bolometric Peak flux 
(FPR) of 1366 BATSE LGRBs (red dots) superposed on the pre- 
dictions of the LGRB world model (black dots) for BATSE LGRB 
detection efficiency. The background grey dots represent the model 
predictions for the entire LGRB population (detected and unde- 
tected). Both BATSE LGRB data and the LGRB world model 
indicate the same (Pearson's) correlation strength of Pfpr-t 90 ~ 
0.67 indicating that the observed joint distribution is not a byprod- 
uct of selection effects in the detection process. 

ence and duration, discarding the traditional definition 
of peak flux as the peak photon counts in 1-second time 
interval. 

4. SUMMARY & CONCLUDING REMARKS 

The primary goal of the presented analysis was to 
model and constrain the luminosity function, spectral 
correlations and energetics of Long-duration class of 
GRBs by exploiting the wealth of information that has 
been buried and untouched in BATSE GRB catalog to 
this date. Below is a summary of steps taken to construct 
the LGRB world model, detailed in Section 

• A sample of 1366 LGRBs is carefully selected by 
the use of Fuzzy C-means clustering algorithm 
based on the spectral parameters and inspection of 
the individual ligh tcur ves of BATSE catalog GRBs 
(Figure [lj Section |2~T 1 . 



• It is proposed that the BATSE LGRB data might 
be very well consistent with being drawn from 
a multivariate log-normal population of LGRBs 
in four comoving-frame LGRB spectral variables: 
the bolometric isotropic 1-second peak luminos- 
ity (Li SO ), the bolometric isotropic emission (Ei SO ), 
the spectral peak energy {E VtZ ) and the duration 
(Tgo z ). Therefore, the observed joint distribution 
of the four LGRB spectral variables: the bolomet- 
ric 1-sec peak flux (Pboi), the bolometric fluence 
(Sboi), the observed spectral peak energy (E p ) and 
the observed duration (Tgo), results from the con- 
volution of the comoving-frame multivariate log- 
normal population with the cosmic rate (i.e. the 
redshift distribution) of LGRBs, truncated by the 
complex LGRB trigger threshold of BATSE Large 
Area Detectors (LADs), as illustrated in Section 
(|2.2|), Equations (p~|45]) . A prescription for mod- 
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eling BATSE LAD detection efficiency is given in 
Appendix [B] 

The LGRB model (Equation [j} is fit to BATSE 
data by maximizi ng t he likelihood function of 
the model (Section |2.3| & Appendix [Cj Equation 
C2). In order to derive the best-fit parameters of 



the model and their corresponding uncertainties, 
an Adaptive Metropolis-Hastings Markov Chain 
Monte Carlo (AMH-MCMC) algorithm is set up 
to efficiently sample from the 16-dimensional like- 
lihood function. The best-fit parameters are ob- 



tained for three LGR B cosmic ra tes: SF R of Hop 



kins fc Beacom| ( |2006[ ), SFR of |Li| fl2008|) fe the pre- 
dieted LURB redshiit distribution ot IButler et all 
(20101 which is consistent with LGRB rate tracing 
cosmic metallicity with a cutoff Z/Z Q ~ 0.2 — 0.5. 

• To ensure the model provides adequate fit to ob- 
servational data, mult ivar iate goodness-of-fit tests 
are presented (Section 2.4 & Figure sets[2j[3j[4j[5]). 



Summarized below, are the principal conclusions 
drawn from the analysis based on the proposed LGRB 
world model: 

1. Energetics: It is expected that the peak bright- 
ness distribution of LGRBs has effective range of 
log(Pb i[erg/s/cm 2 ]) G [—7.11 ± 2.66] correspond- 
ing to P bo[ [erg/s/cm 2 } G [1.70 x 10~ 10 ,3.58 x 
10 -5 ]. This translates to a dynamic 3a range 
- in the comoving- frame - of \og{L iso [erg / s\) G 
[51.53±1.99] corresponding to L iso [erg/s] G [3.46x 
10 49 ,3.38 x 10 53 ]. In addition, a turnover is pre- 
dicted in the differential log (AT) — log(P) dia- 
gram of LGRBs at -P50-300 ~ 0.1 [photons / s / cm 2 } 
in BATSE nominal detection energy range (50 — 
300 [KeV). This is consistent with and further 
extends the apparent flattening in the cumulative 
log(AT) - log(P) diagram of SWIFT LGRBs re- 
ported recently by Butler et al. (20101. 

As for the bolometric fluence and the total 
isotropic emission distributions, a range of 
log(Sboi[erg/cm 2 ]) € [—6.16 ±3.01] corresponding 
to S bo i[erg/cm 2 } G [6.82 x 10- 10 ,7.01 x 10~ 4 ] 
is indicated. This translates to an average dy- 
namic 3cr range - in the comoving-frame - of 
log(E iso [erg]) G [51.93 ± 2.71] corresponding to 



Eyjerg] G [1.66 x 10 49 ,4.46 x 10 54 ] (Sections 3.2 
fcjju] Table [TJ. 

2. Durations & Spectral Peak Energies: The 
comoving-frame spectral peak energies (E PiZ ) 
of LGRBs, is likely well described by a log- 
normal distribution with an average 3a range 
of \og(E PtZ [KeV\) G [2.48 ± 1.12] correspond- 
ing to E p \ z [KeV] G [23,4006] with peak LGRB 
rate at E p ^ z ~ 300 [IfeV]. This translates to 
an effective observer-frame peak energy range of 
log(E p [KeV]) G [1.93 ± 1.22] corresponding to 
E p [KeV] G [5,1427] with peak LGRB rate at 
E p ~ 85 [KeV]. It is also observed that the 
observer-frame Tqq durations of LGRBs peaks at 



T 90 - 30 [s] with a 3a range of T 90 [s] G [1.4,620]. 
This translates to an average 3a range of comoving- 
frame log(Tg 0iZ [s]) G [0.92 ± 1.24] corresponding to 

10 Is 



s] G [0.47, 145] with a peak rate at T 90 z 
(Section pU) Table [TJ. 



3. Spectral Correlations: All four LGRB spectral 
variables: L iso , E iso , E pz & T go . z appear to 
be either moderately or strongly positively corre- 
lated with each other. In particular, a relatively 
strong and ' broad 1 (Pearson's correlation coefficient 
PE iso -E PiZ = 0.58 ± 0.04) and highly significant 
(> 14<r) correlation is predicted between E iso & 
E p _ z of Long-duration class of GRBs. Surprisingly, 
Tgo, z appears to evolve with Li SO & Ei SO such that 
brighter burst s ge nerally tend to have longer dura- 
tions (Section 3.4 Table [I]). This prediction of the 
model, together with the previously reported nega- 
tive correlation of the brightness a nd the duration 
of individual pulses in LGRBs (e.g., Fenimore et al 



1995 Nemiroff 2000) might possibly indicate that 
intrinsically brighter LGRBs contain higher num- 
bers of pulses. 

There is a slight chance that a small fraction (< 50) 
of BATSE LGRBs were misclassified as SGRBs 
by the automated pattern recognition methods ex- 
ploited in this analysis (c.f. Figure 3l center right 
panel). If true, it will most likely affect (if signifi- 
cant at all) the constraints derived on the luminos- 
ity function of LGRBs and the correlation of E PtZ 
with L isn . 



4. Redshift Distribution: The lack redshift informa- 
tion for the BATSE GRBs strongly limits the pre- 
diction power of the presented analysis for the cos- 
mic rate of LGRBs. Nevertheless, based on the 
Markov Chain sampling of the likelihood function 
for the three LGRB redshift distributions consid- 

Figure [9]) , it is observed 
jx not necessarily, 
favors an LGRB rate consistent with cosmic metal- 
licity evolution wi th a cutoff Z/Z® ~ 0.2 — 0.5 (c.f. 



ered here (Section 2.2 
that BATSE data potentially, bu 



Butler et al. 2010 ), with no luminosity-redshift evo- 
lution. 

Assuming LGRBs track Star Formation Rate, only 
a tiny fraction (-2-3) of 1366 BATSE LGRBs 
are expected to have originated from high redshifts 
(z > 5). In the case of an LGRB rate tracing cos- 
mic metallicity evolution (e.g., Butler et al.||2010 ), 
the fraction increases by one order of magnitude to 
~ 2%, corresponding to — 27 bursts out of 1366 
BATSE LGRBs. For comparison, the expected 
fraction of SWIFT fe EXIST LGRBs w ith z £ 5 
are ~ 6% & ~ 7% ( |Butler et aL|[20To| . The dis- 
crepancy is well explained by the fact that both 
SWIFT and EXIST are more sensitive to long-soft 
bursts - characteristic of high redshift LGRBs - 
due to their lower gamma-ray trigger energy win- 
dow ! _comp_ared__to^ATSE_Ijar^ 
(c.f. | Gehrels et al. | 2004||Band et al.pO^ [GTuidlay 



fc Team||2009| ). 



Although fitting is performed for the comoving-frame 
spectral variables, it is notable that the overall shape of 



LGRB World Model 



15 



the resulting observer-frame distribution of the spectral 
variables also resembles a multivariate log-normal (c.f. 
Figures [2] [3j II) |>J [6| 0. In other words, the redshift 
convolution ofthc comoving-frame population distribu- 
tion approximately acts as a linear transformation from 
the comoving frame of LGRBs to the observer frame. 
This is primarily due to the narrow redshift distribution 
of LGRBs - as compared to the width of the LGRBs 
comoving-frame spectral distributions - with almost 90% 
of the popul ation originating from intermediate redshifts, 



rates of these events convolved with complex detection 
thresholds of different instruments used for observations. 



G [1,4.3]. Balazs et al. (2003) provide an elegant dis 



cussion on the potential effects of redshift convolution on 
the observed distribution of LGRBs spectral parameters. 

As implied by the model, there is no evidence for 
a significant population of bright-hard LGRBs that 
could have been missed in BATSE catalog of GRBs. 
Conversely, a large population of low-luminosity with 
moderate-to-low spectral peak energies seem to have 
gone undetected by BATSE LAD detectors. Weather the 
X-Ray Flas hes (XRF), X-ray rich and the sub- luminous 



GRBs (e.g.,[St rohmayc r et al.||1998| |Kippen et al.||2003 
Tikhomirova, et al.|200g| " can be incorporated into a um 



tied class ot events described by a single model, remains 
an open question in this work. At present, these events 
can be either considered as a separate class of cosmo- 
logical events or as the soft-dim tail of the LGRB world 
model that have been mostly out of BATSE detection 
range and missed (c.f. Figure (3) of Kippen et al. (20031 
for a comparison with the predictions ot the LGRB world 
model here in Figures [2jM|2|&J5J). A definite answer 
to this question requires knowledge of the true rate of 
sub-luminous bursts and XRFs based on the observed 



Butler et al. (2010) (B10) have presented an elaborate 
multivariate analysis of SWIFT LGRBs. While provid- 
ing reasonable fit to SWIFT data, the model of B10 is 
primarily aimed at the discovery of the potential sub- 
luminous events that mostly go undetected by gamma- 
ray detectors. Such model, capable of accounting for a 
large population of undetected sub-luminous bursts, is 
proposed at the cost of throwing away parts of informa- 
tion stored in the s pectral parameter s of LGRBs in the 
analysis of B10 (c.f. |Ghirlanda et al.|[2012] ). 

The proposed multivariate log-normal model while re- 
quiring minimal free parameters compared to any other 
statistical model considered in the GRB literature to 
date, provides an accurate comprehensive description of 
the largest catalog of Long-duration GRBs, serving as a 
powerful probe to explore the population properties of 
a large fraction of LGRBs that are missed in spectral 
analyses due to low-quality data or the lack of measured 
redshift, or simply go undetected due to the instrument's 
gamma-ray detector threshold. It is expected that a bet- 
ter definition of peak luminosity that is not limited to 
a specific timescale in the observer-frame of the bursts 
such as th e effective luminosity definition of |Butler| 
et al. (2010) - would result in an improvement in the 



correlations of the time-integrated spectral peak energy 
(E PjZ ) and the isotropic emission (Ei SO ) with the lumi- 
nosity variable (Li so ). Data from future gamma-ray ex- 
periments will enable us to further confirm, improve or 
invalidate predictions of the presented model. 



APPENDIX 
A. GRB CLASSIFICATION 

It is well known that the traditional classification of GRBs based on a sharp cutoff in the observed duration (Tqq) 
distribution of GRBs - usually set at Tgg = 2 [s] is insufficient and can be misleading close to the sharply defined border. 
The apparent long-soft-bright t o short-hard-dim trend observed in the spectral properties of BATSE GRBs (e.g. Figure 
[I] top panel; also Figure (8)]of Shahmoradi & Nemiroff (2010)) necessitates the use of a rigorous classification scheme 
Dased on all available spectral properties ot GRBs, m addition to duration. 

Despite a rich literature on the classification methodologies for Gamma- Ray Bursts (e.g., Hakkil a~et al.|2000a 2003 



and references therein) the choice of a classification method to separate BATSE catalog of GRBs into two subgroups o 
Long & Short -duration with minimal misclassifications remains a difficult task. This is primarily due to the significant 
overlap (or similarity) in all (or some) spectral properties of the two classes of GRBs, in addition to heterogeneity of 
objective functions that might differ considerably from one classification algorithm to another. 

Here, fuzzy (soft) clustering algorithms arc preferred over hard clustering methods, since they provide a probability of 
the event belonging to each specific subgroup, in contrast to hard classifications that return only binary probabilities of 
cither or 1. The ch oice o f fuzzy algorithm greatly facilitates identification of bursts that might have been potentially 
misclassified (Section 



2.1) 



Investigation of different fuzzy algorith ms available in the literature leads us to t wo pr om inent candidat es: The 



SAND method of |Rousseeuw et al.| ( |1996| & the fuzzy C-means discussed by |Dunn| ( [1973] ) fc |Bezdek| ( 1981[ ). While 
fuzzy C-means is specially useful tor cases where subgroups are known to be approximately symmetric, the SAND 
algorithm is superior to C-means for its lack of sensitivity to different subgroup sizes, orientations and asymmetries. 
Nevertheless, the presence of a handful of Soft-Gamma Repeaters (SGRs) in the BATSE catalog - with spectral 
properties comparable to that of LGRBs - results in relatively poor classification by the SAND method as compared 
to C-means. Besides the choice of algorithm, the spectral variables - by which the classification is done - are selected 
such that the resul ting relative sizes of the two SGRB and LGRB populations correspond to those found 



& Nemiroff] ( |2010[ ) through a different approach that the authors believe to be less prone to biases (c.f. Figure |( 13) | 
Table j(4~j] in Ehahmoradi & Nemiroff 2010). 



B. BATSE TRIGGER EFFICIENCY 



Before the LGRB world model of Equation (HI) is fit to BATSE observational data, it is necessary to convolve the 
model with BATSE trigger threshold (Equation^]) . The study of BATSE detection efficiency is well documented in a 
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Fig. 9. — Left: A comparison of the model-predicted BATSE trigger efficiency of 1366 BATSE LGRBs with the nominal trigger efficiency 
estimates of BATSE 4B Catalog. The discrepancy between the two curves is primarily due to different methodologies and GRB samples 
used to derive BATSE Catalog efficiency curvtf ^ and the trigger efficiency of the LGRB world model. Right: The nor mali zed sampling 
distribution of the likelihood function (Equation |C2| from Markov chain for three LGRB redshift distributions (Section |2.2| Equation [JJ . 
The fact that the redshift distribution of B10 generally results in larger likelihoods as compared to the two other, can provide evidence - 
but does not necessitate in any way - that BATSE data favors a LGRB rate tracing cosmic metallicity evolution over SFR. A firm decision 
can only be made by the complete knowledge of BATSE LGRB redshifts. 

a http : //www .batse .msf c .nasa. gov/batse/grb/catalog/4b/tables/4br_grossc . trig_sen 

series of articles 



fc Nemiroff| ( |2011 | 
LGRBs have duration of T90 



by BATSE team (e.g., [Pendleton et al.|1995|p98)|Paciesas et al.|1999||Hakkila|2003] c.f. 
1 1 for further discussion and references). Here, based on the observation that almost all 1 
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references) 

> l[s], the primary trigger time-scale for BATSE LGRBs is assumed to be 1024 [ms]. 
This eliminates the relatively complex dependence of the detection probability (77 in Equation [5]) on the duration of 
the events. The probability of detection for a LGRB is then modeled by the Cumulative Density Function (CDF) of 
log-normal distribution, 



T](detection\lJ,thresh, ^thresh) = 



1 1 f log (P [L iso , E p<z ,z)) - (Ithresh 



-erf 



V2cr 



I 



(Bl) 



where P(Li SO , E p>z , z) is 1-second peak photon flux in the BATSE nominal detection energy range: 50 — 300 [XeF], 
and [ithresh & ^thresh are the detection threshold parameters to be determined by the model. The link between the 
1-sec peak photon flux (P) and the LGRB comoving- fram e variables (Lj sn , E PtZ ) and redshift (z) is provided by fitting 
a smoothly broken power-law known as the Band model ( |Band et al.|[l9 93) to LGRBs differential photon spectra, 



E a 
E? 



exp 



( _ (i+z)(2+ a )E/ Ep ^ E < [E Pt ,/(l + z)] [{a - /3)/(2 + a)] 



otherwise. 



such that, 



P^^iso: Ep Z: z) 



AttDt 



K ' J0.1 



20000/ 



/(!+*) 



E<PdE 



(B2) 



(B3) 



It has been shown by B10 that fixing the the high-energy and low-energy photon indices of the Band model (Equation 

B2) to the corresponding population average a — —1.1, /3 = —2.3 produces only a negligible error of < 0.05 dex in 
the resulting flux estimates. Given the uncertainties in the BATSE LGRB variables , in particular E p estimates, such 
an assumption provides reasonable approximation for the calculation of the peak fluxes. A more accurate treatment, 
however, would be to include possible weak cor relations that are observed bet we en the Band model photon ind ices 
and the spectral peak energies of the bursts (c.f. Shahmoradi & Nemiroff (2010) & Shahmoradi & Nemiroff (2011 1 for 
a discussion of the correlations and potential origins). 

The goodness of the log-normal CDF assumption for BATSE detection efficiency can be checked by a comparison of 
the resulting model predictions with BATSE 1366 LGRBs's distribution of peak fluxes (Figure [2]). The best-fit model 
prediction of BATSE trigger efficiency for Long-duration class of gamma-ray bursts is comparedro the nominal trigger 
efficiency of BATSE 4B catalog for the class of soft-long bursts in Figure ([9j Left Panel). Although the difference 
between the two curves is significant, it does not necessarily imply a contradiction, given the fact that different 
methodologies and GRB samples were used to derive the two efficiency curvetj^J 



1 |BATSE 4B Exposure vs Peak Flux Table 
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C. LIKELIHOOD FUNCTION 



To obtain the joint posterior for the unknown parameters of the LGRB world model of Equation (HI) given BATSE 
data, the likelihood functio n of the model must be, in princ iple, constructed by correctly accounting for uncertainties 
in observational data (e.g., Eddington 1913 Jeffreys 1938). In addition, it is known tha t astronomical survey s at 



low Signal-to-Noise Ratios (SJNR) clos e to survey thres hold can be potentially biased (e.g., Hogg & Turner 1998). A 
Bayesian multilevel methodology (e.g., Hobson et al.|2010 ) can incorporate the above corrections required to construct 
the likelihood function: Under the assumption of normality for the uncertainties of LGRB spectral variables in BATSE 
catalog, each LGRB event - denoted by Oi, standing for the ith LGRB Observation - has the likelihood £, of having 
the true spectral parameters G\ = [Pboi,i, Sboi,i, Ep,i^90,i\ that is described by a 4-dimensional Gaussian Probability 
Density Function (pdf), 



(CI) 



where the location (fJ-ifi) an d the scale (S,,o) parameters of the pdf are to be determined by the model and individual 
events data from BATSE GRB Catalog. Given £j and the LGRB world model (Equation |l]) convolved with BATSE 
trigger efficiency (Equation [5| , the full Poisson likelihood function can be written as, 



£(Data|Model Parameters) = A N exp -A / K obs {0) dO 

\ J Ospace 

N r 

X JJ / Tlcosmicifiifi, £j,o, ^|Model Parameters)£i(Oi|/4o, Si,o) dO, (C2) 

i—l J Ospace 

where N = 1366, and A is a factor that properly normalizes the cosmic and the observed rates (7Z C osmic & TZ-obs)- 
In the absence of the prior knowledge (i.e. lZ C osmic, as in the case here), the Empirical Bayes approach can provide 
an alternative solution, in which, an ad hoc estimate of the model parameters based on the observed data serves 
as the prior for same data at the second level of analysis. Calculation of the normalization factor A, involves an 
integration of the LGRB world model over the 5-dimensional space of LGRB spec tral variables and redshift, with a 
complex integration limit set by BATSE trigger efficiency modeled in Equation Bl as a function of Li SO , E v%z & z. In 
addition, since almost no redshift information is available for BATSE catalog of GRBs, the probability for each LGRB 
has to be marginalized over redshift, for which a range of z € [0.1, oo] is considered. These integrations make the 
maximization of the likelihood for all unknown parameters an extremely difficult task, perhaps challenging the current 
computational technologies. Moreover, it has been known that GRB fluences and durations are likely und erestimated 



Hakkila et aLl|2000b[ |2003|). Such bias 



close to detection threshold due to the so-called fluence-duration bias (e.g 
makes the overestimation-correction of the fluence and duration as prescribed by 
before the fluence-duration bias effects are well quantified. The algorithms for calcu 

to result in less biased measurements - even down to detector threshold - with negligible uncertainties 



Hogg & Turner| |"l998 ) unjustified 
ating peak fluxes, however, a ppear 



et al 



B van; 



(.e.g., 



Stern 



Given the above lines of reasoning and the computational limitations, the u ncertainties on three BATSE" 



LGRB variables: Pboh Sboi, T90 are excluded from the likelihood function (Equation C2). The uncertainties on the E j 



estimates of BATSE LGRBs are, however, significant compared to the three former variables and must be incorporated 
in the calculations of the likelihood. Nevertheless, it was realised aft er likelihood maximization that the exclusion of 
the E p uncertainties - by fixing E p values to the 'bisector' estimates of Shahmoradi & Ncmiroff (2010) - results in only 
negligible (<C lc) changes in the mo del's best-fit para meters. In general, the use ot the bisector line of Ordinary Least 



Squares (OLS) regression lines (e.g., Isobe et al.|1990[) for estimation purposes is unfavored due to lack of a maximum 



likelihood interpretation. The special case of BATSE LGRBs here, however, turns out to be an exception 



In addition to cosmological time-dilatio n correction 



i t is common practice to make an energy correction to the 
temporal variables of SWIFT GRBs (e.g., Gehrels et al. 2006 Butler et al.||2010 ), such as T go , when translating the 
variable from the observer-frame to the comoving-lrame. For B Al SE LGRBs, this energy-correction is likely negligible, 
given the fact that the T90 durations are calculated based on the total photon counts in the BATSE LAD energy range 
.|1993[ Fishman et al.|19"94 1 : 20 — 2000 KeV which can be practically considered as bolometric. 
;cted that such an energy correctic 



(e.g., Kouveliotou et al. 



Nevertheless, it is expected that such an energy correction, if needed, would slightly relax the tight correlation of the 
comoving-frame duration {Tg 0z ) with the total isotropic emission (E iso ) and the peak isotropic luminosity (L iso ). 

The join t po sterior distribution of the model parameters are obtained by the maximizing the likelihood function of 
Equation ( C2 ) convolved with non-inf ormative unifo rm prior on the location parameters and the standard choice of 



Jeffreys prior on the scale parameters ( | Jeffreys]|1946[ ) . In order to efficiency sample from the 16-dimensional posterior 
density function, a variant of Marko v Chain Monte Ca rlo (MCMC) methods known as Adaptive Metropolis-Hastings 
(AMH-MCMC) is employed (e.g., |Haario et~aO2001[ ). The choice of an Adaptive (vs. classical) MH algorithm is 



very important, since the model parameters exhibit strong covariance (Table p5|). To reduce the simulation runtime, 
all algorithms including AMH-MCMC are implemented in Fortran, which isby far the fastest and most efficient 
programming language for many intensive scientific calculations and number crunchingr] In addition, the numerical 
integration in the definition of the luminosity distance (Equation[3j) - encountered on the order of J> 10 8 times during full 

3 The Ideal HPC Programming Language 
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MCMC sampling - is greatly simplified by the analytical approximation method of Wickramasinghe & Ukwatta (2010). 
Due to the intrinsic sequential character of MCMC sampling methods, the parallelization of simulation algorithms tor 
a single Markov chain is impractical or at best inefficient. Nevertheless, to increase the MCMC sample size and more 
importantly, to ensure convergence to the global (vs. local) extremum, the chain is initiated simultaneously at 10 
random starting points in the parameter space on a 12-corcs desktop CPU. In general, convergence and good mixing 
occurs within the first few thousands of iterations (burn-in period) given a suitable initial guess for the covariance 
matrix of the proposal distribution - here chosen to be multivariate Gaussian. The resulting mean and la standard 
deviations of the model parameters are tabulated in Table ([!]). 

D. LGRB MONTE CARLO UNIVERSE 

The prediction power and consistency of the presented LGRB world model - based on BATSE data - can be easily 
checked against observational data from current and future gamma-ray experiments, in particular Fermi satellite. All 
it takes is to construct a Monte Carlo universe of LGRBs, based on the best fit parameters of the LGRB world model in 
Table and compare the outcome with observational data. Although straightforward, the steps for such simulation 
and comparison are summarized bellow: 

1. A random redshift for each simulated LGRB is drawn from the redshift distribution of Equation Q with 
parameters taken from Table ([I]). It is recommended to repeatedly randomly draw the set of model parameters 
from the full Markov Chain sample^] instead of fixing the parameters to the mean values reported in Table ([I]). 

2. The four LGRB spectral variables: Li SO , Ei SO , E PtZ , Tgo tZ are randomly drawn from a 4-dimensional log-normal 
distribution with location (p) and scale (i.e., the covariance matrix: E) parameters constructed from fitting 
results in Table 0. This can be easily and quickly done by noting that a multivariate log-normal distribution 
is equivalent to a multivariate Gaussian distribution in the logarithmic space of the above spectral variables, 



logO= I log(L iso ),log(E iso ),log(E PtZ ),log(T 90tZ )j, (Dl) 

such that the 4-dimensional log-normal density function , CAf, of Equation ([I]) can be exactly replaced by a 
4-dimcnsional Gaussian distribution, 

Af(logO\(l,z) = (2n)- k / 2 \j:\- 1 / 2 x exp ( - (logO- fl'JT 1 (log O-fifi), (D2) 



for which, the cosmic LGRB differential rate (TZ cosm i C ) of Equation ([!]) will be, 

Ucosmic = (D3) 

dlog{L lso ) dlog(E iso ) dlog(E PtZ ) dlag(T 90iZ ) dz 1 

3. The above Monte Car lo universe of LGRBs can be then measured according to the instrument's detection 
efficiency (Section |2.2| Equation [5]). For the case of BATSE LAD detectors, the trigger efficiency can be 
modeled according to the prescription in Appendix |B| 



The author owes special thanks to Robert J. Nemiroff (NASA GSFC Astrophysicist & professor of Physics at 
Michigan Technological University) for his generous support and help through numerous discussions over more than 
a year, from the early stages of developing the basic ideas behind this work to the core of the analysis and final steps 
of manuscript preparation. 

I would like to thank William H. Press (professor of Computer Science and Integrative Biology at The University 
of Texas at Austin) and Lars Koesterke (at Texas Advanced Computing Center) for helpful comments on some 
statistical and computational aspects of this work. I am also very grateful to Swadesh Mahajan (professor at Institute 
for Fusion Studies) & Richard Hazeltine (professor and chair at the Department of Physics, The University of Texas 
at Austin) for their generous support in the final stages of preparing this manuscript. 

This work would have not been accomplished without the vast time and effort spent by many scientists and engineers 
who designed, built and launched the Compton Gamma-Ray Observatory and were involved in the analysis of GRB 
data from BATSE Large Area Detectors. 



Available at https : //sites . google . com/site/amshportal/ — r—, j- — — 7^ — c 73 3~n 

I £ — — S — S : . £ U research/aca/in-the-news/lgrb-world-model 



LGRB World Model 

TABLE 2 1366 BATSE catalog triggers classified as LGRBs 



19 



L11 gg cl 




mggCl 






r 1 Vi <T£TCH* 


i rigger 


i rigger 


' I Vi (rn'/^Y" 
lrl SS cr 


T 1"1 (TITt^'V 

lrl 66 cr 


' 1 Vi crcyc* v 


r 1 Vi cfcr&v 


' 1 Vi rprfiiv 

i rigger 




105 


107 


109 


110 


111 


114 


121 


130 


133 


143 


148 


160 


171 


179 


204 


211 


214 


219 


222 


223 


226 


228 


235 


237 


249 


257 


288 


332 


351 


394 


398 


401 


404 


408 


414 


451 


465 


467 


469 


472 


473 


493 


501 


516 


526 


540 


543 


548 


549 


559 


563 


577 


591 


594 


606 


630 


647 


658 


659 


660 


673 


676 


678 


680 


685 


686 


690 


692 


704 


717 


741 


752 


753 


755 


761 


764 


773 


795 


803 


815 


816 


820 


824 


825 


829 


840 


841 


869 


907 


914 


927 


938 


946 


973 


999 


1009 


1036 


1039 


1042 


1046 


1085 


1086 


1087 


1114 


1120 


1122 


1123 


1125 


1126 


1141 


1145 


1148 


1150 


1152 


1153 


1156 


1157 


1159 


1167 


1190 


1192 


1196 


1197 


1200 


1204 


1213 


1218 


1221 


1235 


1244 


1279 


1288 


1291 


1298 


1303 


1306 


1318 


1382 


1384 


1385 


1390 


1396 


1406 


1416 


1419 


1425 


1432 


1439 


1440 


1446 


1447 


1449 


1452 


1456 


1458 


1467 


1468 


1472 


1492 


1515 


1533 


1540 


1541 


1551 


1552 


1558 


1559 


1561 


1567 


1574 


1578 


1579 


1580 


1586 


1590 


1601 


1604 


1606 


1609 


1611 


1614 


1623 


1625 


1626 


1628 


1642 


1646 


1651 


1652 


1653 


1655 


1656 


1657 


1660 


1661 


1663 


1664 


1667 


1676 


1687 


1693 


1700 


1701 


1704 


1711 


1712 


1714 


1717 


1730 


1731 


1733 


1734 


1740 


1742 


1806 


1807 


1815 


1819 


1830 


1883 


1885 


1886 


1922 


1924 


1956 


1967 


1974 


1982 


1989 


1993 


1997 


2018 


2019 


2035 


2047 


2053 


2061 


2067 


2069 


2070 


2074 


2077 


2079 


2080 


2081 


2083 


2087 


2090 


2093 


2101 


2102 


2105 


2106 


2110 


2111 


2112 


2114 


2119 


2122 


2123 


2129 


2133 


2138 


2140 


2143 


2148 


2149 


2151 


2152 


2156 


2181 


2187 


2188 


2189 


2190 


2191 


2193 


2197 


2202 


2203 


2204 


2207 


2211 


2213 


2219 


2228 


2230 


2232 


2233 


2240 


2244 


2252 


2253 


2254 


2267 


2276 


2277 


2287 


2298 


2304 


2306 


2309 


2310 


2311 


2315 


2316 


2321 


2324 


2325 


2328 


2329 


2340 


2344 


2345 


2346 


2347 


2349 


2362 


2367 


2371 


2373 


2375 


2380 


2381 


2383 


2385 


2387 


2391 


2392 


2393 


2394 


2405 


2419 


2423 


2428 


2429 


2430 


2432 


2435 


2436 


2437 


2438 


2440 


2441 


2442 


2443 


2446 


2447 


2450 


2451 


2452 


2453 


2458 


2460 


2472 


2476 


2477 


2482 


2484 


2495 


2496 


2500 


2505 


2508 


2510 


2511 


2515 


2519 


2522 


2528 


2530 


2533 


2537 


2541 


2551 


2560 


2569 


2570 


2581 


2586 


2589 


2593 


2600 


2603 


2606 


2608 


2610 


2611 


2619 


2620 


2628 


2634 


2636 


2640 


2641 


2660 


2662 


2663 


2664 


2665 


2671 


2677 


2681 


2688 


2691 


2695 


2696 


2697 


2700 


2703 


2706 


2709 


2711 


2719 


2725 


2727 


2736 


2749 


2750 


2751 


2753 


2767 


2770 


2774 


2775 


2780 


2790 


2793 


2797 


2798 


2812 


2815 


2825 


2830 


2831 


2843 


2848 


2850 


2852 


2853 


2855 


2856 


2857 


2862 


2863 


2864 


2877 


2880 


2889 


2890 


2891 


2897 


2898 


2900 


2901 


2913 


2916 


2917 


2919 


2922 


2924 


2925 


2927 


2929 


2931 


2932 


2944 


2945 


2947 


2948 


2950 


2951 


2953 


2958 


2961 


2980 


2984 


2985 


2986 


2990 


2992 


2993 


2994 


2996 


2998 


3001 


3003 


3005 


3011 


3012 


3015 


3017 


3026 


3028 


3029 


3032 


3035 


3040 


3042 


3055 


3056 


3057 


3067 


3068 


3070 


3071 


3072 


3074 


3075 


3076 


3080 


3084 


3085 


3088 


3091 


3093 


3096 


3100 


3101 


3102 


3103 


3105 


3109 


3110 


3115 


3119 


3120 


3127 


3128 


3129 


3130 


3131 


3132 


3134 


3135 


3136 


3138 


3139 


3141 


3142 


3143 


3153 


3156 


3159 


3166 


3167 


3168 


3171 


3174 


3177 


3178 


3193 


3212 


3217 


3220 


3227 


3229 


3237 


3238 


3241 


3242 


3245 


3246 


3247 


3255 


3256 


3257 


3259 


3267 


3269 


3276 


3279 


3283 


3284 


3287 


3290 


3292 


3301 


3306 


3307 


3319 


3320 


3321 


3322 


3324 


3330 


3336 


3339 


3345 


3347 


3350 


3351 


3352 


3356 


3358 


3364 


3369 


3370 


3378 


3403 


3405 


3407 


3408 


3415 


3416 


3436 


3439 


3448 


3458 


3465 


3471 


3472 


3480 


3481 


3485 


3486 


3488 


3489 


3491 


3493 


3503 


3505 


3509 


3511 


3512 


3514 


3515 


3516 


3523 


3527 


3528 


3552 


3567 


3569 


3588 


3593 


3598 


3608 


3618 


3634 


3637 


3648 


3649 


3654 


3655 


3658 


3662 


3663 


3664 


3671 


3717 


3733 


3740 


3745 


3765 


3766 


3768 


3771 


3773 


3776 


3779 


3788 


3792 


3800 


3801 


3805 


3807 


3811 


3814 


3815 


3819 


3840 


3843 


3853 


3860 


3869 


3870 


3871 


3875 


3879 


3886 


3890 


3891 


3892 


3893 


3899 


3900 


3901 


3903 


3905 


3906 


3908 


3909 


3912 


3913 


3914 


3916 


3917 


3918 


3924 


3926 


3929 


3930 


3935 


3941 


3954 


4039 


4048 


4095 


4146 


4157 


4216 


4251 


4312 


4350 


4368 


4388 


4556 


4569 


4653 


4701 


4710 


4745 


4814 


4939 


4959 


5080 


5255 


5304 


5305 


5379 


5387 


5389 


5407 


5409 


5411 


5412 


5415 


5416 


5417 


5419 


5420 


5421 


5423 


5428 


5429 


5433 


5434 


5447 


5450 


5451 


5454 


5463 


5464 


5465 


5466 


5470 


5472 


5473 


5474 


5475 


5476 


5477 


5478 


5479 


5480 


5482 


5483 


5484 


5486 


5487 


5489 


5490 


5492 


5493 


5494 


5495 


5497 


5503 


5504 


5507 


5508 


5510 


5512 


5513 


5515 


5516 


5517 


5518 


5523 


5524 


5526 


5530 


5531 


5538 


5539 


5540 


5541 


5542 


5545 


5548 


5551 


5554 


5555 


5559 


5563 


5565 


5566 


5567 


5569 


5571 


5572 


5573 


5574 


5575 


5581 


5585 


5589 


5590 


5591 


5593 


5594 


5597 


5601 


5603 


5604 


5605 


5606 


5608 


5610 


5612 


5614 


5615 


5617 


5618 


5621 


5622 


5624 


5626 


5627 


5628 


5632 


5635 


5637 


5640 


5644 


5645 


5646 


5648 


5654 


5655 


5667 


5697 


5704 


5706 


5713 


5715 


5716 


5718 


5719 


5721 


5723 


5725 


5726 


5729 


5731 


5736 


5773 


5867 


5890 


5955 


5983 


5989 


5995 


6004 


6082 


6083 


6090 


6098 


6100 


6101 


6102 


6103 


6104 


6111 


6113 


6115 


6118 


6119 


6124 


6127 


6128 


6131 


6137 


6139 


6141 


6147 


6151 


6152 


6154 


6158 


6159 


6165 


6167 


6168 


6176 


6186 


6188 


6189 


6190 



Continued on Next Page 



20 



Shahmoradi 
TABLE 2 - Continued 



Frigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigger 


Trigge: 


6194 


6198 


6206 


6222 


6223 


6225 


6226 


6227 


6228 


6233 


6234 


6241 


6242 


6243 


6244 


6249 


6266 


6267 


6269 


6270 


6271 


6272 


6273 


6274 


6279 


6280 


6283 


6285 


6288 


6295 


6298 


6300 


6303 


6304 


6305 


6306 


6308 


6309 


6315 


6317 


6319 


6320 


6321 


6322 


6323 


6328 


6329 


6330 


6334 


6335 


6337 


6339 


6344 


6345 


6346 


6349 


6351 


6353 


6355 


6369 


6370 


6375 


6380 


6388 


6390 


6395 


6396 


6397 


6399 


6400 


6404 


6405 


6408 


6409 


6413 


6414 


6419 


6422 


6425 


6435 


6437 


6440 


6444 


6446 


6448 


6450 


6451 


6453 


6454 


6472 


6487 


6489 


6490 


6498 


6504 


6519 


6520 


6521 


6522 


6523 


6525 


6528 


6529 


6531 


6533 


6534 


6536 


6538 


6539 


6544 


6546 


6550 


6551 


6552 


6554 


6557 


6560 


6564 


6566 


6576 


6577 


6578 


6582 


6583 


6585 


6587 


6589 


6590 


6592 


6593 


6598 


6600 


6601 


6602 


6605 


6610 


6611 


6613 


6615 


6616 


6619 


6620 


6621 


6622 


6625 


6629 


6630 


6631 


6632 


6642 


6648 


6649 


6655 


6657 


6658 


6665 


6666 


6670 


6672 


6673 


6674 


6676 


6678 


6683 


6686 


6694 


6695 


6698 


6702 


6707 


6708 


6720 


6745 


6762 


6763 


6764 


6767 


6774 


6782 


6796 


6802 


6814 


6816 


6830 


6831 


6853 


6877 


6880 


6882 


6884 


6891 


6892 


6903 


6911 


6914 


6917 


6930 


6935 


6938 


6963 


6987 


6989 


7000 


7012 


7028 


7030 


7064 


7087 


7108 


7110 


7113 


7116 


7130 


7147 


7164 


7167 


7170 


7172 


7178 


7183 


7185 


7191 


7206 


7207 


7209 


7213 


7219 


7228 


7230 


7247 


7250 


7255 


7263 


7285 


7293 


7295 


7298 


7301 


7310 


7318 


7319 


7322 


7323 


7328 


7335 


7343 


7357 


7358 


7360 


7369 


7371 


7374 


7376 


7377 


7379 


7381 


7386 


7387 


7390 


7403 


7404 


7429 


7432 


7433 


7446 


7451 


7452 


7457 


7460 


7464 


7469 


7475 


7477 


7481 


7485 


7486 


7487 


7488 


7491 


7493 


7494 


7497 


7500 


7502 


7503 


7504 


7509 


7515 


7517 


7518 


7520 


7523 


7527 


7528 


7529 


7532 


7533 


7535 


7548 


7549 


7550 


7551 


7552 


7560 


7563 


7564 


7566 


7567 


7568 


7573 


7575 


7576 


7579 


7580 


7587 


7588 


7597 


7598 


7603 


7604 


7605 


7606 


7607 


7608 


7609 


7614 


7615 


7617 


7619 


7625 


7630 


7635 


7638 


7642 


7645 


7648 


7654 


7656 


7657 


7660 


7662 


7677 


7678 


7683 


7684 


7688 


7695 


7701 


7703 


7705 


7707 


7711 


7727 


7729 


7741 


7744 


7749 


7750 


7752 


7762 


7766 


7769 


7770 


7780 


7781 


7785 


7786 


7788 


7790 


7794 


7795 


7798 


7802 


7803 


7810 


7818 


7822 


7825 


7831 


7835 


7838 


7840 


7841 


7843 


7845 


7858 


7862 


7868 


7872 


7884 


7885 


7886 


7888 


7900 


7902 


7903 


7906 


7918 


7923 


7924 


7929 


7932 


7934 


7936 


7938 


7942 


7948 


7954 


7963 


7968 


7969 


7973 


7976 


7984 


7987 


7989 


7992 


7994 


7997 


7998 


8001 


8004 


8008 


8009 


8012 


8019 


8022 


8026 


8030 


8036 


8039 


8045 


8049 


8050 


8054 


8059 


8061 


8062 


8063 


8064 


8066 


8073 


8075 


8084 


8086 


8087 


8098 


8099 


8101 


8102 


8105 


8110 


8111 


8112 


8116 


8121 















Notes: Spectral data for these triggers are available in BATSE 4B Sz Current Catalogs The spectral peak energy (E p ) estimates of the above 
triggers and the rest of 2130 BATSE Catalog GRBs are provided by [Shahmoradi &z Nemiroff (2010|. The full conditional E p probability density 
functions are available for download at https://sites.google.com/site/amshportal/research/aca/iri-the-riews/lgrb-world-model 
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TABLE 3 

Correlation Matrix of the parameters of the LGRB World Model, 

FOR THE MEDIAN CASE OF AN LGRB COSMIC RATE TRACING STAR FORMATION RATE OF Li (2008) 
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